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ABSTRACT 

The Jeans equations give the second moments or stresses required to support a stel- 
lar population against the gravity field. A general solution of the Jeans equations for 
arbitrary axisymmetric scale-free densities in flattened scale-free potentials is given. 
A two-parameter subset of the solution for the second moments for the self-consistent 
density of the power-law models, which have exactly spheroidal equipotentials, is ex- 
amined in detail. In the spherical limit, the potential of these models reduces to that 
of the singular power-law spheres. We build the physical three-integral distribution 
functions that correspond to the flattened stellar components. 

Next, we attack the problem of finding distribution functions associated with the 
Jeans solutions in flattened scale-free potentials. The third or partial integral intro- 
duced by de Zeeuw, Evans & Schwarzschild for Binney's model is generalised to thin 
and near-thin orbits moving in arbitrary axisymmetric scale-free potentials. The par- 
tial integral is a modification of the total angular momentum. For the self-consistent 
power-law models, we show how this enables the construction of simple three-integral 
distribution functions. The connexion between these approximate distribution func- 
tions and the Jeans solutions is discussed in some detail. 
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1 1 INTRODUCTION 

Evidence both from hydrodynamical modelling of ellipti- 
cals (e.g., Binney, Davies & Illingworth 1990; van der Marel 
1991) and the kinematics of the stellar halo of the Milky 
Way (e.g., Norris 1986; Morrison, Flynn & Freeman 1990) 
indicates that the phase space distribution functions (DFs) 
of elliptical galaxies and the haloes of spiral galaxies depend 
on three isolating integrals of motion. 

General flattened potentials support fewer than three 
global isolating integrals. In triaxial systems, only the en- 
ergy E is globally conserved. In axisymmetric potentials, 
the angular momentum component parallel to the symme- 
try axis L z is also an invariant. Nonetheless, most stars do 
admit a third integral of motion that is a generalisation of 
the total angular momentum L (c.f., Saaf 1968; Innanen & 
Papp 1977; Lupton & Gunn 1987). The existence of chaotic 
trajectories prevents the extension of an analytic third in- 
tegral to all the orbits. So, building three-integral DFs is 
often difficult in axisymmetric potentials. The problem is 
still more severe for triaxial systems (Schwarzschild 1979). 

The natural shape of ellipticals and the haloes of spiral 
galaxies may well be triaxial. Nonetheless, it is still both 
useful and realistic to consider the somewhat simpler prob- 
lem of how to construct three-integral axisymmetric stellar 
systems. Recently, de Zeeuw, Evans & Schwarzschild (1996, 
hereafter ZES) showed how to find a third, partial, integral 
of good accuracy for thin and near-thin tubes for oblate 



scale-free potentials with flat rotation curves. They used the 
partial integral to build three-integral DFs and identified a 
class of physical solutions to the Jeans equations. The aim of 
this paper is to extend their methods and results to arbitrary 
axisymmetric scale-free potentials, with particular emphasis 
on the models with spheroidal equipotentials known as the 
power- law models (Evans 1994, hereafter E94). 



The paper is arranged as follows. In Section 2, the Jeans 
approach is applied to yield the complete solution for the 
stresses or second order velocity moments that can support 
flattened populations of stars in scale-free potentials. A par- 
ticularly simple two-parameter solution is identified for the 
self-consistent power-law models. Section 3 investigates the 
spherical limit. New families of three-integral DFs are pre- 
sented and used to isolate the Jeans solutions of physical in- 
terest. With the insight gained from this simpler problem, we 
begin the harder job of constructing three-integral DFs for 
flattened scale-free potentials in Section 4. The third or par- 
tial integral introduced by ZES for the scale-free logarithmic 
potential (e.g., Binney & Tremaine 1987; see also Richstone 
1980) is generalised to arbitrary axisymmetric scale-free po- 
tentials. We show by detailed numerical orbit integrations 
that it is a good integral for the thin and near-thin tube or- 
bits, although it is not so well-conserved for fat tube orbits. 
Three integral DFs are then built in Section 5 and shown to 
correspond to a subset of the earlier Jeans solutions. 
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2 2 THE JEANS APPROACH 



2.2 2.2 General solution of the Jeans equations 



In this section, we use standard spherical coordinates 
(r,9,<j>), with 9 measured from the axis of symmetry and 
<f> the azimuthal angle. 



2.1 2.1 Scale-free potentials and densities 

Scale-free axisymmetric potentials have the form 



$ = _- 



(2.1) 



where g(9) is an arbitrary function that describes the shape 
of the equipotentials. The circular velocity v c i IC in the equa- 
torial plane varies like r _/3//2 . So, models with /3 < have ris- 
ing rotation curves, whereas models with /3 > have falling 
rotation curves. When (3 = 0, the scale-free power-law po- 
tential becomes logarithmic, and we recover the special case 
already studied by ZES. The relevant range of (3 is from 1 
corresponding to the outer Keplerian envelopes of star clus- 
ters to —2 corresponding to the shallow cusps of boxy ellip- 
ticals. Henceforth, we use units in which the velocity scale 
fo is unity. 

It is often useful to consider tracer populations of stars 
moving in an external gravity field rather than the self- 
consistent stellar density generated by Poisson's equation. 
So, let us take the density to have the general scale-free 
form: 



P = 



h(9) 



rV +/3/2 (#)' 



(2.2) 



Here, 7 is a constant, which prescribes the radial fall-off of 
the density, while h(9) is an arbitrary function. Both g(9) 
and h(9) arc symmetric with respect to the equatorial plane, 
i.e., g(n-9) = g(9) and h(n-9) = h(6). When p is the self- 
consistent density, then 7 = 2 + (3 and h(8) and g(9) are 
related by (4ttG = 1) 



h(9) =(l-p)g 2 (0) + \g(0)g'(0) cot 9 

i S W(0)(l + i/3). 



(2.3) 



When 7 7^ 2 + /3, the density (2.2) is greater than the density 
associated with the potential (2.1) either at large radii (7 < 
2 + 0) or at small radii (7 > 2 + 0). However, we are often 
interested in the dynamics in particular regimes - such as 
the outer reaches or the cusp - and so this is not a grave 
drawback, provided the difficulty occurs outside the regime 
under scrutiny. 

One pleasing choice for the arbitrary function g(6) is 



9(0) 



cos 2 e 



(2.4) 



so that the equipotentials are similar concentric spheroids 
with axis ratio q. These are recognised as the power-law 
models introduced in E94. The associated self-consistent 
density is of the form (2.2), with 

h(0) = Q{(l-/3g 2 )sin 2 + [2-Q(l + /3)]cos 2 6»}, (2.5) 



with Q = q ' 
h{9) = 1-/3. 



In the spherical limit, q = Q = 1 and so 



The potential (2.1) and density (2.2) have the felicitous at- 
tribute of scale-freeness. Their properties at radius r' = kr 
are just a magnification of those at radius r. Scale-free den- 
sity distributions may have DFs that are not scale-free. For 
example, scale-free spheres can possess DFs built according 
to the instructions provided by Osipkov (1979) and Merritt 
(1985). These have an anisotropy radius, at which the prop- 
erties of the velocity dispersion tensor change. However, let 
us assume that the associated DFs are also scale-free (e.g., 
Richstone 1980; White 1985; ZES), then the stresses have 
the following form: 



p(v r ) = 



P{ve) 



Fi(9) 



r f3 + ~fg2 + 

F 3 (0) 

r fj+-, n 2+f3 



(9)' 



, 1 F 2 (9) 
P^ve) = r/3+7 „ 2+/3 



F 4 (9) 



(2.6) 



r f3+ ln 2+f3 



where Fi, F2, F3, and F4 are functions of 9, fulfilling the 
following conditions 



Fi(tt- 
F 2 (tt- 



Fi(9), 
-F 2 (9), 



F A (it- 
F 4 (tt- 



F 3 (9), 
F 4 (9). 



(2.7) 



As discussed in ZES, these constraints guarantee that the 
stresses are symmetric with respect to the equatorial plane. 
Substitution of the Ansatz (2.6) into the Jeans equations 
(ZES, equation 2.4) reduces them to two coupled first-order 
ordinary differential equations: 



F2(9) + 



cot 9 



(2 + /%'(#) 



9(0) 



F 2 (9) 



(2.8a) 



+ (2-/3- 7)Fi(0) - F S (9) - F 4 (9) = -h(9), 



Fs(9) + 



cot 9 



(2 + /%'(#) 



9(0) 



Fz{9) 



+ (3-/3- -y)F 2 (9) - F 4 (9) cot0 = - 



h(9)g'(9) 
MO) ' 



(2.8b) 



The four functions Fi(9), . . . , F 4 (9) are therefore subject to 
the above two restrictions. We are at liberty to pick two of 
the functions arbitrarily and solve (2.8) for the other two. 
We choose to prescribe Fi(9) and F 2 (9). Proceeding as in 
ZES, we use equation (2.8a) to eliminate F 4 (9) from (2.8b). 
It is straightforward to integrate the resulting first order 
differential equation for J3 and then to substitute the result 
in (2.8a) to obtain F 4 . We find: 



F 3 (9) = 1(9) + J(9) + F 2 (9) cot 9, 



F 4 (9) = h(9) - 1(9) - J {ff) 

-(P + 1 -2)F 1 (9)+g f3+2 (9)^ § 



F 2 (9) 



9 p+2 (9) 



(2.9a) 



(2.9b) 



with 



sin 2 J gP+ 2 (0) ] 




2g(0) } 



Simple Three- Integral Scale-Free Galaxy Models 3 



and 



,13+2 



(0) 



[F 2 (9)-F 1 (9) cote] 



(2.11) 



Hence, the stresses can be found by evaluation of the inte- 
grals for 1(9) and J (6). It follows that 

F 3 (0) = F 4 (0) = ift(O) - i(/3 + 7 - 2)ft(0). (2.12) 

This means that (v%) = (v 2 ,) on the minor axis, as is required 
by elementary symmetry arguments (e.g., Bacon 1985). Let 
us note that the case 7 + (3 = 2 is special, as then the 
integral J (9) drops out of the solution (2.9). The only self- 
consistent model (7 = 2 + /3) for which this happens has (3 — 
and 7 = 2. This completes our derivation of the general 
solution of the Jeans equations for galaxies with scale-free 
potentials (2.1) and scale-free densities of the form (2.2). 
The condition that the principal components of the stress 
tensor are positive definite is given in section 2.4 of ZES. 

2.3 2.3 The two-integral limit: / = f(E, L 2 ) 

Let us consider briefly the special case corresponding to a 
two-integral DF / = f(E,L 2 ), in which the stellar veloci- 
ties have no preferred direction in the meridional plane, i.e., 
(v 2 ) = (vg) and (v r vg) = 0. It is straightforward to find the 
associated stresses solution by taking 



Fi(9) = F z (9), 



F 2 (9) = 0. 



(2.13) 



Substitution of (2.13) in equations (2.9a) and (2.11) results 
in a Volterra integral equation for F\_ which is solved to give 



Fi(9) 



(0) 



dd [cot 9 - 



9 0+2 (6) 



9(0) ■ 

MO)- 



h(9). (2.14) 



Table 1. Coefficients for the Jeans solution (2.20), 
with Q = q~ 2 . 

A 1 Hi(l + /3-/3q 2 ) + H 2 (l + /3) + ±(l-/3) 

Ci QH 2 (1 + (3) - (3Hi + ±Q[2 - Q(l + /3)] 

B 2 #i(l + /3) 

A 3 1(1 -ft- fa 2 Hi 

C 3 H 1 + ±Q(2-Q(l + /3)) 

A 4 Q(i + /3) - i(i + 3/3) -ff 2 (l + /3)(l + 2/3) 
+f3Hi{q 2 (l + 2(3)- 2(1 + f3)} 

C A Qi? 2 (l + /3)-/3i?i + iQ[2-Q(l + /3)] 



Here, 1(9) + J(9) is 



1(9) + J(9) 



1 1 



2 1 + 



P ■ 2, 

— sm ( 



1+/3 



H ig 2 g(9) 



+ 



1Q[2- 



— — cos 9, 



(2.18) 



1 + 13 



while Hi and H 2 are constants. Of course, this is a very re- 
stricted subset of the general solution (2.9) - instead of two 
free functions Fi(9) and F 2 (9), we have merely two free pa- 
rameters Hi and H 2 . They have a simple physical interpre- 
tation. Hi = means that the velocity ellipsoids are aligned 
with the cylindrical coordinate system, while H 2 = means 
that they are aligned with the spherical coordinate system. 

Substitution of the choice (2.17) in the general solution 
(2.9) yields 



F z (6) = 1(9) + J(9) + H 2 g(9) cos 2 9, 
F A (9) = h(9) - 1(9) - J (9) - 2f3Fi(9) 

+ H 2 [Q cos 2 9 -sin 2 9 - 2/3(1 - Q) sin 2 9 cos 2 9]. 



(2.19) 



Substitution of this expression in equation (2.9b) provides 
F 4 (9) as: 



F 4 (9) = h(9) + (l-p -j)Fi (9). 



(2.15) 



These results can be checked using Hunter's (1977) solution 
of the Jeans equations for general f(E,L 2 z ) models. 

2.4 2.4 A two-parameter solution for the 
self-consistent power-law models 

Let us now specialise to the specific case of the self-consistent 
power-law models. We take g(9) and h(9) as given in (2.4) 
and (2.5), while 7 = 2 + (3. Evans & de Zeeuw (1994) showed 
that when the stresses have the simple form 



P(Vj) = 



aiR 2 + bjRz + CiZ 2 



(2.16) 



(R 2 + z 2 /q 2 ) 2 +P ' 

then the line-of-sight projected second moments - such as 
the dispersions in the radial velocities and the proper mo- 
tions - can be explicitly evaluated as elementary functions. 
As a matter of choice, we elect to study the subset of the 
general solution (2.9) that possesses this fetching property. 
This implies choosing 



Fi (9) = 1(9) + J (9) + Hi + H 2 g(9) sin 2 9, 
F 2 (9) = H 2 g(9)sm9cos9. 



(2.17) 



This seems complicated, but when we write the Jeans solu- 
tion out explicitly in cylindrical coordinates, it is simply: 



p{vr) 



p(vrv z ) 



P{Vz) 



P(vl) 



AiR 2 + Cn 



(1 + 0)(R 2 + z 2 /q 2 Y+f i '' 

B 2 Rz 

(1 + P)(R 2 + z 2 /q 2 ) 2 +f } '' 

A 3 R 2 + C 3 z 2 
(l+(3)(R 2 + z 2 /q 2 ) 2 +l 3 ' 

A4.R 2 + CiZ 2 
(l + P)(R 2 + z 2 /q 2 ) 2 +f' 



(2.20) 



where the relationships between the coefficients Aj , Bj , Cj 
and our two anisotropy parameters Hi and H 2 are given in 
Table 1. Lastly, let us remark that the two-parameter Jeans 
solution (2.20) admits an extension to cored models that is 
discussed in Appendix A. 



2.5 2.5 Limits on the two-parameter solution 

The principal stresses p(v 2 ) and p(v 2 ) of the two-parameter 
Jeans solution for the power-law models are 



P(vl) 



F-(0) 



rP + fn 2 



g 2 {oy 



P{vl) 



F+(0) 



rP+ia 2 



g 2 {ey 



(2.21) 
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Figure 1. The (Hi , /i2)-planc of two-paramctcr anisotropic so- 
lutions of the Jeans equations for the self-consistent power-law 
models with /3 = —0.25 and q = 0.87. The hatched area corre- 
sponds to the solutions of the Jeans equations with non-negative 
stresses p{v^), p(v z ), and p{v^). The requirement that the prin- 
cipal components of the stress tensor p(v^) and p(v 2 ) be positive 
definite leads to two further small areas in the bottom left quad- 
rant being deleted. The line Hi = gives the cylindrically aligned 
models, while the line H2 = gives spherically aligned models. 



(2.22) 



with 

F±{8) = 1(9) + J '(9) + \Hx + \H 2 g(9) 

± ^H 2 - 2HiH 2 g(6) cos 26» + H 2 g 2 (6). 

The tilt angle Q of the velocity ellipsoid in the meridional 
plane is the misalignment with respect to the spherical polar 
coordinate surfaces. It is: 

tan2e - 2 Il - H 2 g(6) sin26 

tan 20- ^ - Hl _ H2g{0)cos29 - V-W 

Necessary (but not sufficient) requirements for positive 
stresses are that p(v%), p(v 2 ) and p(u|) are non-negative. 
This implies the following conditions on Hi and H2: 

-\Q[2-Q(l + l3)]<Hi, 

±Q(l+f3)-l<(l + f3)H 2 -f3q 2 Hi, 

-i(l-/?)<(l + /J)ff a 

+ (l+(3-f3q 2 )Hi, 

Q(l + /3)- i (1 + 3/3) > H 2 (l+/3)(l + 2(3) 

-/3[q 2 (l + 2/3)-2-2/3]Hi. 

These define a hatched region in the (Hi, H2) -plane illus- 
trated in Figs. 1 and 2 for models with rising and falling 
rotation curves respectively. Note that when /3 < 0, the 
hatched region extends to infinity along the H 2 = axis. 
When f3 > 0, there exists a limiting value of Hi, namely 
Q(l — 2/3 + q 2 )/(4/3), beyond which the Jeans solutions be- 
come negative. The case f3 = discussed in ZES is interme- 
diate, in the sense that the upper and lower boundaries of 



(2.24) 



Figure 2. As Fig. 1, but for fi = 0.18 and q = 0.9. The hatched 
area now closes, and it is not possible to find Jeans solutions with 
arbitrarily large radial anisotropy. Note that the largest value 
of Hi consistent with positive stresses is Q(l — 2/3 + <jr 2 )/(4/3). 
This occurs when the lines Ci = and A4 = intersect. 
Adding in the additional condition that the solution is spher- 
ically aligned (H 2 = 0) gives the more stringent condition 
pHi <Q(l-iQ(l + /?)). 

the hatched region are exactly horizontal. 

Necessary and sufficient conditions for positive stresses 
are that the principal components p(vx) and p{v 2 ) (or, 
cquivalently, their sum and product) are non-negative. This 
problem can be solved using the same methods as in section 
2.7 of ZES. The additional regions that must be discarded 
from Figs. 1 and 2 are those simultaneously satisfying the 
four inequalities written out in detail in Appendix B. In 
fact, this results in the removal of only two small areas near 
the two corners of the hatched region that lie in the third 
quadrant of Figs. 1 and 2, respectively. They lie away from 
the spherically aligned Jeans solutions (H2 = 0), for which 
we will construct approximate DFs in Section 5. As these 
additional, tiny forbidden areas are not particularly impor- 
tant for the purposes of this paper, we have not complicated 
Figs. 1 and 2 by marking them. 

The position of the two-integral solution (2.13) is 
marked in Figs. 1 and 2 - it lies at the origin in the (Hi,H2)- 
plane. As shown in E94, it is generated by a positive definite 
DF provided the flattening satisfies q 2 > i(l + /3). Note that 
this constraint can also be deduced from the first equation of 
(2.24) on putting Hi — 0. It is clear from Figs. 1 and 2 that 
the set of Jeans solutions for the self-consistent power-law 
models is large. The rest of the paper addresses the ques- 
tion: which solutions are physically relevant and correspond 
to positive DFs? 



3 3 SPHERICAL POTENTIALS 

Now let us suppose the potential is exactly spherical. The 
beauty of this assumption is that it enables the construction 
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of flattened components with simple three-integral DFs and 
triaxial kinematics. This is because the total angular mo- 
mentum L is an exact integral of motion in a spherical po- 
tential well, in addition to the angular momentum compo- 
nent about the symmetry axis L z . As first realised by White 
(1985), such DFs give realistic descriptions of the kinematics 
of tracer populations of stars in the outer reaches of galaxies. 



3.1 3.1 General Jeans Solutions and Distribution 
Functions 

If the gravity field is spherically symmetric, then g(9) = 1 
and (2.1) reduces to the potential of the singular power-law 
spheres (E94). The work in the previous section allows us to 
deduce the general solution for the stresses associated with 
a flattened density p — r *h(8). It is of the form (2.6) with 
g(6) = 1. The functions Fi(9) and F2(6) are arbitrary, while 
F 3 (6) and F A (Q) follow from (2.9) with g(6) = 1. However, 
in a spherical potential, we must have (v r ve) = (see e.g., 
ZES). So, only the solutions with Fziff) = are physical! 
They have the form 



F 3 (6) = I(6)+J(0), 

F4(0) = fc(0)-/(0)-J(0)-(/9 + 7-2)fi(0), 



(3.1) 



and 



1(6) 



sin 2 9 



h(9) sin 9 cos 9d9, 



(3.2) 



Fi(0) sin 6» cos 9d9. 



Here, F\(9) - or the angular variation of the radial velocity 
dispersion - is arbitrary. Once it and h(9) have been chosen, 
the other second moments are fixed. When /3 + 7 = 3, we 
have (v 2 ) = (v 2 ) + (v 2 ) + (v 2 ,) = r~ f3 . A special case of this 
rule was noted earlier by Maoz & Bekenstein (1990), who 
pointed out that if the velocity dispersion is independent of 
both r and 9, then the potential is logarithmic (/3 = 0) and 
the density profile must fall like r~ 3 (i.e., 7 = 3). 

Now h(9) and Fi(9) are even functions symmetric about 
6 = 7r/2. They may be expanded as a series of even powers 
of sin6l (see e.g., Mathews & Walker 1964; ZES): 



h{6) = ^ a n sin 2n 6, Fi(0) = ^ &„ sin 2n (9. 



(3.3) 



Our aim is to find DFs generating the Jeans solution (3.1) 
corresponding to each pair h(9),F\(9). First, we note that 
the stellar density law r -7 sin 2n 9 can be reproduced by non- 
negative DFs of the form 



f m , n (E,L 2 ,L 2 z ) = Vm , n L 2m Ll n \E\ a , 
where 

(2-p)(m + n) 7 3 
P P 2' 

and m + n > —1, 2n > —1 and 



7 + 2m + 2n > 



0, 

(3{m + n + 



if (3 < 0, 
if > 0. 



(3.4) 



(3.5) 



(3.6) 




Figure 3. The (Hi , i?2)-plane of two-paramctcr anisotropic so- 
lutions of the Jeans equations for the singular self-consistent 
power-law spheres with rising rotation laws (/3 < 0). The hatched 
region indicates the area with everywhere positive stresses. The 
rilled circle corresponds to the isotropic model (Hi = 0, H2 = 0). 
The physical solutions are those that are spherically aligned - 
these lie on the thickened, bold line H2 = 0, subject only to the 
condition that Hi > — i(l — /3). 



Note that we have written m and n so that our notation is 
consistent with White (1985) and ZES, but m and n are not 
necessarily integers. The normalization constant r) m ,n is 

|^|2m//3+2n/ / 3+ 7 / / 3 



0F2 m +"+ 3 / 2 r(m + n + 1)5(1/2, n + 1/2) 



r(-«) 

T(-7//3 - 2m/ (3 - 2n/f3) ' 

r( 7 //3 + 2m/ /3 + 2n/(3 + 1) 
r(a + l) 



if f3 < 0, 



if P > 0. 



(3.7) 



Here, B(x,y) — T(x)T(y) /T(x + y) is the beta function. 
By integrating over all velocity space, the corresponding 
stresses are 



p{v 2 ) = 
P(ve) = 



1 



2 m + 2n + P + 7 



r/3 + 7 



m + n + 1 



P(v r ), 



(3.8) 



(2n + \)p(vl). 



The angular dependence of the density is independent of m 
for fixed 7 and n. The same is true of the stresses. However, 
the anisotropy ratios (v 2 )/(v 2 ) and (v 2 ,) / ' (v 2 ) do vary with 
m. Whenever 7 + P = 3, then (v 2 ) + (v 2 ) + (v 2 /) = f 2 irc . 
In the limit ft — > 0, the components reduce to those stud- 
ied in White (1985), Gerhard (1991) and ZES. The compo- 
nents for P > have previously been discussed by Kulessa 
& Lynden-Bell (1992) and de Bruijne, van der Marel & de 
Zeeuw (1996). Using the method of linear superposition of 
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Figure 4. As Fig. 3, but for the self-consistent power-law spheres 
with falling rotation laws (/3 > 0). Note that the physical solutions 
(corresponding to the thickened, bold line) are more restricted 
than in Fig. 3. They encompass those parts of the line H2 = 0, 
subject to -1(1 - P) < Hi < 1(1 - 13)/ 13. 

the components (3.4), the general solution of the Boltzmann 
equation becomes 



f(E,L 2 ,L 2 z ) 



(3.9) 



Here, the A m ,n are constants specifying the fraction con- 
tributed by each component. To build DFs that recover the 
angular variation h(ff) of the density and Fi(0) of the ra- 
dial velocity dispersion of our Jeans solution, we must insist 
that: 



E 



A, 



E 



2m+2n+ /3 + 7 



(3.10) 



for n = 0, 1, 2, . . .. There are many more unknowns A m . n 
than constraints a n ,b n . So, there are many ways to define 
three-integral DFs that reproduce the required density and 
velocity dispersions. For specific choices of h(9), examples 
of such DFs are given in Kulessa & Lynden-Bell (1992) and 
de Bruijne, van der Marel & de Zeeuw (1996). 

As an aside, we remark that it is possible to generalise 
(3.4) to produce triaxial stellar distributions with triaxial 
kinematics in spherical potentials. In view of its ready ap- 
plication to the kinematics of the stellar halo of the Milky 
Way Galaxy, this development is discussed in more detail in 
Appendix C. 



3.2 3.2 A one-parameter solution for 
self-consistent spherical models 

Let us now consider the limiting behavior of the two- 
parameter Jeans solutions for the self-consistent power-law 
spheres. Not just the potential, but also the density, is now 



spherical. Using h(0) = 1-/3, g(8) 
find from equations (2.17) that 



1 and 7 = 2 + /3, we 



Fi(0) = 



1 - /3 + 2#i 



2(1 +P) 
F 2 {6) = # 2 sin (9 cos (9, 
1 - P - 2/?#i 



+ H 2 sin 2 0, 



F 3 (9) 
F 4 (6) 



2(1 + 0) 
1-P- ipHx 
2(1 + /?) 



+ Hi cos 



+ # 2 [cos 2 - sin 2 0(1 + 2/3)]. 



(3.11) 



These anisotropic stresses for the self-consistent, singular 
power- law spheres are positive in the region in the (Hi,H 2 )- 
plane defined by (cf. eq. [2.24]) 



■\{i-P)<m, 

-i(l-/3)<ff 2 (l + /3) 



pm, 

■5(1-0) < H 2 (l + P) + H u 
1(1-/3) > # 2 (l + /3)(l + 2/3)+/3ffi 



(3.12) 



which is illustrated in Figs. 3 and 4 for power-law spheres 
with rising (/3 < 0) and falling (/3 > 0) rotation laws. As we 
have seen, all physical solutions must have F 2 = because of 
the properties of individual orbits. On requiring that H2 = 
0, we are left with a one-parameter solution 



Fi(0) 



1-P + 2H 1 



2(1 

F 3 (0) = F 4 (0) 



1-P- 2/3#i 
2(1 + p) ' 



(3.13) 



which is indicated as the thick solid line in Figs. 3 and 4. It 
follows from the derivation in Section 3.1 that each of these 
Jeans solutions corresponds to a DF of the general form (3.9) 
with n — 0, so that 



f(E,L 2 ) = ^2A m , i)Vm ,oL 2m 



\E\ 



2(m+l)//3- 



'1/2 



(3.14) 



with rimfi given in equation (3.7). The coefficients are sub- 
ject to the constraints 



J2 Am,o = 1-/3, J2 



A-r, 



1 



/3 + 2 J ffi 



m+P+1 



1 + P 



(3.15) 



Again, many DFs are possible. The simplest is obtained by 
taking only one component in the series (3.14): 



f(E,L 2 
with 



L 2m \E\ 2{r ' 



■.+l)//3-m-l/2 



-2ffiQ3 + l) 
1 + 2i3i - P ' 



(3.16) 



(3.17) 



Of course, Hi is a measure of the anisotropy of the model. 
Hi — —1(1-/3) corresponds to the circular orbit model 
(no radial velocity dispersion) and Hi = to the isotropic 
model. When /3 < 0, then Hi — » 00 and the model ap- 
proaches the radial orbit model. When /3 > 0, then the 
most radially distended model possible (see Fig. 4) has 
Hi = (1 — P)/(2p). This corresponds to the critical value 
of m = — 1, for which (3.4) ceases to exist. The DF (3.16) 
is non-negative and so physical for each of these permitted 
values of Hi. 
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Table 2. Accuracy of the partial integral I3 for an E3 power-law 
model with a falling rotation curve (/3 = 0.18,5 = 0.9). Meridional 
cross sections of the orbits are displayed in Fig. 5. 



q Ro zo 


L z 


^3, max 


To 

I 5,min 


57< 

u 1 s 


L 2 

max 


L 2 

mm 


5L 2 


Thin tubes in 


E3 gal 


ixies 












A 0.9 .5232 .2384 


.541 


.3678 


.3677 


.0001 


.3677 


.3534 


.0143 


B 0.9 .3794 .4268 


.378 


.3671 


.3660 


.0011 


.3660 


.3236 


.0424 


C 0.9 .1889 .5367 


.182 


.3657 


.3630 


.0027 


.3630 


.3004 


.0626 


D 0.9 .0472 .5667 


.045 


.3651 


.3617 


.0034 


.3617 


.2935 


.0682 


Fat tubes in E3 galaxies 












E 0.9 .4800 .4853 


.378 


.3499 


.3208 


.0291 


.3484 


.2889 


.0595 


F 0.9 .6600 .4529 


.378 


.2583 


.2208 


.0375 


.2531 


.2094 


.0437 


G 0.9 .8400 .2425 


.378 


.1666 


.1555 


.0111 


.1650 


.1539 


.0111 



Table 3. Accuracy of the partial integral Iz for an E3 power-law model 
with a rising rotation curve (f3 = —0.25, q = 0.87). 



q Ro zo 




-^3, max 


^3,min 


5/3 


1-2 

max 


L 2 . 

mm 


5L 2 


Thin tubes in 


E3 galaxies 












A 0.87 .5852 .2573 


.541 


.3677 


.3674 


.0003 


.3674 


.3493 


.0181 


B 0.87 .4334 .4677 


.378 


.3653 


.3628 


.0025 


.3627 


.3094 


.0533 


C 0.87 .2220 .6003 


.182 


.3600 


.3527 


.0073 


.3527 


.2753 


.0774 


D 0.87 .0562 .6395 


.045 


.3566 


.3477 


.0089 


.3477 


.2640 


.0837 


Fat tubes in E3 galaxies 












E 0.87 .4800 .4961 


.378 


.3678 


.3438 


.0240 


.3676 


.2955 


.0721 


F 0.87 .6600 .4927 


.378 


.3037 


.2444 


.0593 


.3010 


.2225 


.0785 


G 0.87 .8400 .3423 


.378 


.2079 


.1742 


.0337 


.2027 


.1666 


.0361 



4 4 THE BOLTZMANN APPROACH 

The Boltzmann approach emphasises the primacy of dis- 
tribution functions (DFs). In this section, we investigate 
which of the two-parameter Jeans solutions can correspond 
to three-integral DFs for the flattened power-law potentials. 
To do this, we must first extend the theory of partial inte- 
grals introduced in ZES. 



4.1 4.1 Classical and partial integrals 

The collisionless Boltzmann equation for the DF / can be 
cast into the form 

= v r A + veB, (4.1) 
with 

A 1 9 f ^< 2 ^ 2 d ^^ d f 2d f 2 d f (A K 

A =* r fr+to+ V +- r fr)M- V9 W- V +d%' (4 ' 2) 

and 

,df 9^ df 2 fdf df\ 
B - 2 M-wd^ +v * cote {d^~d^)- (4 ' 3) 

Of course, any axisymmetric potential has two classical in- 
tegrals, the energy E = $ + \{v 2 + vj + v 2 ,) and the z- 
component of the angular momentum L z — rv^sinO. With 
the exception of special potentials possessing additional 
symmetries in phase space, axisymmetric potentials do not 
have exact, globally defined third integrals (e.g., Lynden- 
Bell 1962; de Zeeuw 1985; Evans 1990). This is awkward for 



stellar dynamics, as the observational data require three- 
integral DFs. One way round this problem is suggested in 
ZES. We look for partial integrals which have good accuracy 
for some orbital families, rather than the global integrals 
which we know in general do not exist. ZES constructed 
such a partial integral for the thin and near-thin tubes in a 
scale- free model with a flat rotation curve (Binney's model). 
Here, our aim is to understand how this partial integral ex- 
tends to the entire scale-free family. 

In the spherical limit, the potential has an exact third 
integral, namely the total angular momentum L 2 — r 2 (vg + 



ZES introduced a modification 



L 2 i 3 (9) = r 2 (vl + vl)i 3 (6), 



(4.4) 



and showed that a suitable choice for Binney's model was 
h{0) = sin 2 9 + Qcos 2 0. Inserting our ansatz 73 into the 
Boltzmann equation (4.1), we find that 



A =0, 



B =\r%{9) 



2 . 2 
v e + V4, - 



i',(9)gi+?/ 2 (9) 



The thin or near-thin tubes have the property that 
<9<& 1 



dr rPgP/ 2 {9)' 



(4.5) 



(4.6) 



is fulfilled. This essentially states that anywhere on a thin 
tube, v r is much smaller than the other two components. 
This itself implies that, in the meridional plane, the thin 
tubes lie very nearly on circles. This is certainly a good 
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approximation in moderately flattened models (like E3) . The 
partial integral therefore satisfies the collisionless Boltzmann 
equation if 



1 d«3 1 dg 
h~dB = ~g"d~8' 

Straightforward integration gives 

is(0)=g(0), 

and so the third integral is just 
h = r 2 (vg +vl)g(0). 



(4.7) 



(4.8) 



(4.9) 



This is exactly the same as deduced by ZES for the specific 
case of Binney's model. We emphasise, though, that the only 
assumption involved in deriving (4.8) is that of scale-freeness 
of the potential. The claim is that this is a good partial 
integral for thin and near-thin tubes in any galaxy with a 
scale-free potential. 

4.2 4.2 Numerical integrations 

The accuracy of the partial integral is now investigated by 
adaptive Runge-Kutta techniques for the specific cases of 
the self-consistent scale-free power-law galaxies (i.e., g(0) is 
given by (2.4)). The results of such computations are listed 
in Tables 2 and 3, while Fig. 5 shows cross-sections of the 
orbits in the meridional plane. To ease comparison with the 
results of ZES, the energy surface is always normalised to 



E = 



P 



■ exp 



P 



(4.10) 



This ensures that the limiting angular momentum of the 
circular orbit is e -1 ^ 2 . For the purposes of illustration, 
we concentrate on one model with falling rotation curve 
{P = 0.18, q = 0.9, E = -5.5811) and one model with ris- 
ing rotation curve ((5 = -0.25, q = 0.87, E = 4.0268). The 
rationale behind the choices of the equipotential axis ratio 
q is that the galaxies have the same asymptotic ellipticity 
of E3 (using formula (2.9) of E94). The values of f3 are sug- 
gested by taking the means of the asymptotic logarithmic 
gradients of the rotation curves in the samples of Casertano 
& van Gorkom (1991). In Tables 2 and 3, the orbits are 
identified by their starting values Ro and zq, with v r and 
v$ zero there. The maximum and minimum values of 73 and 
L 2 are recorded, together with their fluctuation. Figure 5 
shows cross-sections of the orbits for (3 = 0.18 (the diagram 
for f3 — —0.25 is very similar). The thin tubes are labelled 
A, B, C, and D. For all these thin tubes, the partial integral 
is well-conserved. The lower section of Tables 2 and 3 shows 
the results for thicker tube orbits. Now the variation in the 
partial integral is larger - but it is still more accurately con- 
served than the total angular momentum. The partial inte- 
gral is not exact, even for the thin tubes. At any angular 
momentum L z , there is a unique infinitesimally thin tube 
with energy (4.10). Fig. 6 shows the fluctuation in 73 (ex- 
pressed as a percentage) along the thin tubes for the three 
values p — 0.18,0.0 and —0.25. This illustrates graphically 
that the fluctuations in the partial integral are smaller for 
the models with falling rotation curves. Note that, even for 
the thin tubes, the partial integral is not exact. 

Can the partial integral accurately distinguish between 
different orbital families or - equivalently - different tori 



q = 0.9 
|S= 0.18 




Figure 5. Cross sections of the orbits of Table 2 with a merid- 
ional plane (R, z) for the scale-free power-law potential q = 0.9 at 
energy E = —5.5811. The zero-velocity curve (for L z = 0) is indi- 
cated by the ZVC. Orbits A - D are thin tubes. Orbits E - G arc 
fat (or thick) tubes that fill the indicated areas. (Cross-sections 
of the orbits in Table 3 look very similar). 




Figure 6. The fluctuation in the partial integral for the infinites- 
imally thin tubes (expressed as a percentage) is plotted against 
the angular momentum for three values of f). Even for the thin 
tubes, the partial integral is not exact, although it is very good. 
The diagram clearly illustrates that the partial integral worsens 
in accuracy as /3 diminishes. (All the models have q = 0.9). 

in phase space? One way to establish this is to investigate 
the phase space structure using Poincare surfaces of section 
(see e.g., Gutzwillcr 1990). Fig. 7 shows surfaces of section 
for the model with a falling rotation curve (/3 = 0.18, q = 
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5 5 THREE-INTEGRAL DISTRIBUTION 
FUNCTIONS 

We now exploit the partial integral for the scale-free power- 
law models to build simple three-integral DFs by superpos- 
ing basic component solutions of the collisionless Boltzmann 
equation. First, we develop some preliminary results con- 
cerning our components in Section 5.1. Then, in the next 
section, we pass to the actual construction of the approxi- 
mate analytical three-integral DFs by superposition of the 
components. 




Figure 7. Poincare surfaces of section for the model with 
P = 0.18 and q = 0.9. The energy hypersurface is normalised to 
E = —5.5811. The upper and lower panels show the the (R, pjj) 
and (R,p z ) cross-sections of the energy hypersurface at the equa- 
torial plane for L z = 0.5 and at L z = 0.1 respectively. Each 
orbit is marked with the fluctuation in the partial integral rep- 
resented as a percentage. At high L z , the orbits arc mainly thin 
and near-thin tubes. As is evident from the upper two panels, 
the partial integral is well-conserved by all the orbits. At low L z , 
minor orbital families trapped around the 1 : 1 and 3 : 4 reso- 
nances (the 'reflected banana' and 'reflected fish' families) make 
their appearance. The lower two panels show the partial integral 
is poorly conserved for these families. Near the thin-tube orbit, 
though, the conservation is good. 



0.9, E — —5.5811). The upper and lower panels show the 
(R,Pr) and (R,p z ) cross-sections at the equator for large 
and small L z respectively. In each case, the orbits are marked 
with the percentage fluctuation of the partial integral. At 
high L z , there is the reassuring picture of good conservation 
of the partial integral all across the surface of section. The 
orbits are all tubes. At low L z , the resonances associated 
with the 'reflected banana' and 'reflected fish' families (see 
e.g. Lees & Schwarzschild 1992; E94) make their appearance. 
The partial integral is poorly conserved on the islands in 
the lower panels of Fig. 7. Only near the thin-tube orbits 
and near the orbits confined to the equatorial plane is the 
preservation of the partial integral at all good in the lower- 
panels. It is useful to contrast this behaviour with that of 
the classical integral L 2 , which does not distinguish between 
any of the orbits on the lower panel of Fig. 7. We conclude 
then that I3 is a useful discriminant between orbits of the 
same E and L z . 

Although we have only presented numerical evidence for 
the power-law models, which have spheroidal equipotentials, 
nonetheless the derivation given in Section 4.1 suggests that 
our partial integral is good for thin and near-thin tubes in all 
axisymmetric scale-free potentials. It would be interesting to 
substantiate this claim for other scale-free galaxy models - 
such as those studied by Toomre (1982) for instance. 



5.1 5.1 Components 

The two exact integrals I\ = \E\, I2 = L z and the par- 
tial integral is can be used to construct component DFs 
f(Ii, hjh)- Let us note that all such DFs correspond to 
spherically aligned Jeans solutions, i.e., 



(v r Vg) = (vgVtf,) = (v r V<j,) = 0. 



(5.1) 



This general result follows directly from the forms of the 
integrals. In particular, it depends on the absence of ve- 
locity cross-terms in the partial integral. Is it possible to 
improve the accuracy of the partial integral by incorporat- 
ing velocity cross-terms, which generate misalignment of the 
velocity ellipsoid with the spherical polar coordinates? We 
have explored this question with numerical orbit integra- 
tion. It appears that some slight improvement of the partial 
integral is possible, but only at the cost of increasing math- 
ematical complexity of the partial integral. We believe that 
the coefficient of any cross-term must be small, and so any 
misalignment of the stresses from spherical polars is also 
limited. 

To build three-integral DFs, we shall find it useful to 
study the properties of two sets of components. The first 
set is just the product of powers of the integrals (cf., Frickc 
1951; White 1985; ZES) 



fk,: 



jk jn rm 
1 1 1 2 1 3 ■ 



The density generated by (5.2) is 

r (2-/3)(m+n)-/3(fc + 3/2) gin 2n 



pk 



'[9(0)7 



(5.2) 



(5.3) 



[c/(#)P (fe+m+ ™ )/2+3/3/4 
To reproduce the r-dependence of the density, we must have 
, (2-/3)(m + n) , 2 1 



/3 (5 2' 

so that (5.3) simplifies to 
S n . m sin 2 " e 



Pn,'i 



with 



r 2 +/3 [ 5 (6»)] 1+n+ ' 3 / 2 ' 



(5.4) 



(5.5) 



Sm.n — 



0F2 m+n+3/2 r(m + n + 1)5(1/2, n + 1/2) 

\fi\\ + 2(m + n + l) / fj 

r(-l-2(m + n + l)//3) 



r(l/2-2//3-(2-/3)(m + n)//3)' if P < °' (5.6) 



T(l/2 + 2//3 + (2 - f3)(m + n)/0) 
T(2 + 2(m + n + l)//3) 

The second set of components is 

gk,n,m = h l2(h + Ch) m ■ 



if f3 > 0. 



(5.7) 
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Here, C is a constant, which we assume to be small and check 
a posteriori. Again, to reproduce the r-dependence of the 
density, k is restricted to (5.4). The density corresponding 
to (5.7) is 

- 2na r ,m{2n + 1) sin 2 (6»)- 



Pn,m — 



Sn.m 

r 2 ^ [<7(60] 1+n +' V2 



1 + C- 



2(n + l) g{9) 



(5.8) 



This follows on taking the Taylor expansion of (5.7) and then 
using (5.5) twice. In the spherical limit, we shall show that 
C = 0. Therefore, both sets of components (5.2) and (5.7) 
are axisymmetric generalisations of the components (3.4) in 
the spherical limit. 

5.2 5.2 Three-Integral DFs 

DFs that can reproduce the power-law density are given by 
restricting attention to just the n — and n = 1 components 
in (5.5) and (5.8). We choose to investigate DFs of the form: 



/ — y ' Ao, m gk,0,m + j4.1,m/fc,l,r, 



7 2// 3+ (2- /3)m //3-l/2 (j3+c/2) , 



(5.9) 



+ A 1 , m / 1 2/ ^ +(2 -' 3)(m+1)/ ^ 1/2 7 2 73 m 



where the Ao <m and A\ 
generates the density 

1 



are unknown amplitudes. This 



P = 



r 2+/3 5 2+/3/2(#) 



g{6) ^ A^ m So, r , 



+ sin 2 9 ( ^mCA , m S , m + Ai^Si^ 



(5.10) 



If the DF only depends on globally defined classical inte- 
grals, then it obviously satisfies the Jeans equations. This is 
not guaranteed for DFs depending on any partial or approx- 
imate integrals. So, we must also compare the kinematics 
of the DFs (5.9) with the solutions of the Jeans equation 
derived earlier in Section 2. By straightforward integration 
over velocity space, we discover that the DFs (5.9) generate 
the second moments: 

1 [~1 /^\ X Ao.mS() : rn 



p(vl) 



p{vi)-- 



r 2 + 2/3 ff 2 + /3( 6 ,) [2^ ' (l + m + p) 

m 



•E 



(2- 



(3) 



1 



(1 + m)A 0:m S 0: , 

r 2 + 2/3 ff 2+/3( 6 ,) [2^"' ^ (l +m + f3 ) 
m 

, m(l + m)A 0}rn So,m 



(l + m + /3) 



(2 + m)Ai, m Si, 



+ (2 + m + /3) 

m 

1 (1+™)^™% 



(5.11) 



r 2+2f) g 2+fiy2^"> ^ (l + m + /3) 
m 

+ ±(7 sin 9y — y —- - — — 

m 

(2 + m)A hm S 1 



+ 1 sin 2 



'E 



(2 + m + P) 



As expected, this is a spherically aligned solution and all the 
cross-terms of the velocity dispersion tensor vanish. 

Now, the density of the scale-free power-law models can 
be expanded as: 



r 2+/3 5 2+/3/2(#) 



(2-Q(l + /?)) ff (0) 



+ (/3 + 2)(Q-l)sin 2 6» 



(5.12) 



The spherically aligned Jeans solution can be deduced from 
(2.20) by putting H 2 = 0. We obtain: 



P(Vr) 



p(v 2 e) 



P(vl) 



(1 + /3)- 1 



r 2+2/3 ff 2+/3(0) 



\{2-Q{l + p))+H iq 2 ]g{9) 



+ (Q-l)(l + /3)(i + Hig 2 )sin 2 6» 



(1 + /3)" 1 



r 2+2f3 g 2+l3(Q} 



[i(2-Q(l + /3))-/3g 2 tf 1 ] 5 (0) 



+ ±(Q-l)(l + /3) sin 2 



(5.13) 



(1 + /3)" 1 



r 2+2/3 ff 2+/3(0) 



[i(2-Q(l + /3))-/3Q 2 Jf 1 ] 5 (e) 



+ (Q-l)(l + /3)(| -2/3g 2 iJi)sin 2 6i 



We must now choose the amplitudes Ao ;m and Ai >m and 
the constant C so that the densities (5.10) and (5.12) and 
the kinematics (5.11) and (5.13) coincide. The calculation is 
straightforward and reveals that 

C=-|(Q-1). (5.14) 

This verifies that C is small, at least for moderate Battenings 
This justifies the first-order Taylor expansions used in our 
analysis. For example, if the flattening is E3, then q ~ 0.9, 
and so the next term in the Taylor expansion is 0(C 2 ) which 
is ~ 8%. The constraints concerning the unknown ampli- 
tudes Ao, m and Ai ;Tn are most compactly written by intro- 
ducing the weights 



A , m S , m = (2 - Q(l + P))W , m , 
Ai, m Si, m = (/? + 2)(Q-l)Wi, m . 
Then, there are three conditions on the Wo, m , namely 



(5.15) 



\ - W , m 1 2q 2 H 1 

2^i + m + p 1 + /3 + (l + /3)(2-Q(l + /3))' ( 5 - 16 ) 



^ mWo,: 



2/3g 2 #i 



3[2-Q(l + /?)]■ 

/ft 

There are two conditions on the Wi fm , namely 

E^^ + ^by^- 



E 



Wi,, 



1 + 2g 2 J ff 1 



(5.17) 



2 + m + P 2 + (3 



Provided we choose at least three of the weights Wo, m and 
at least two of the weights Wi, m , these linear equations 
(5.16) and (5.17) can always be solved. Of course, it is not 
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guaranteed that the DFs formed by superposing the com- 
ponents are positive definite. This is a hard matter to check 
analytically, but easy to do on a finite grid with the com- 
puter. In practice, we did not find positivity to be a real 
constraint here. We have verified numerically that some of 
the DFs for the model with (3 = 0.18, q = 0.9 are indeed 
positive definite. In particular, we examined the four cases 
Hi — ±0.1 and ±0.3, with the non-vanishing weights cho- 
sen as Wo,o,Wo,io,Wo,2o,Wi ; o and Wi,2o- It seems that, at 
least for moderate Battenings and anisotropies, many of the 
three-integral DFs constructed according to this recipe are 
everywhere positive in phase space. 

The DFs (5.9) therefore seem to show that the spher- 
ically aligned Jeans solutions (H2 = 0), which are marked 
by the bold lines in Figs. 1 and 2, are physical. This con- 
clusion certainly holds good for tangentially anisotropic so- 
lutions with Hi < 0. When Hi becomes large and positive, 
the stress tensor is radially distended. This is caused by the 
presence of eccentric radial orbits in the model, for which 
our partial integral is not a good invariant. In this limit, 
it is unwise to put undue faith in our DF. In the spherical 
limit, C = Wi t m = Ai >m — 0. This means that the second 
sum in (5.9) vanishes. The partial integral 73 reduces to L 2 
and so we recover the DFs given in Section 3 of this paper, 
as we should. 

There is one further, somewhat special, limit that de- 
serves mention. When (3 = 2, the scale-free power-law po- 
tential is not useful for modelling galaxies, as it generates 
negative density through Poisson's equation. Nonetheless, it 
is perfectly reasonable to consider the properties of orbits in 
this supplied force field. The case is exceptional because the 
Hamilton- Jacobi equation separates (e.g., Landau & Lifshitz 
1960; Lynden-Bell 1962). So, in this limit, there is an exact, 
global third integral of the form 



'3,g = 



2/ 2 . 



1 

W) 



(5.18) 



Now, our partial integral (4.9) does not reduce to this global 
integral when [3 = 2. This, however, is no problem. Any 
thin tube, which is a one-dimensional manifold, must pos- 
sess conserved quantities additional to the globally defined 
integrals. The partial integral valid for thin tubes (4.9) need 
not therefore reduce to (5.18). When [3 = 2, we see that 



I 3 = 2/ 3 ,gS(6>) + 1. 



(5.19) 



Any change in the partial integral A/3 therefore satisfies 



AI 3 = 2/3, g A 5 (0), 



(5.20) 



We see that this change must be small for the thin tubes, 
because /3 ig w from (4.6). Note, too, that it follows from 
the form of (5.18) that the velocity ellipsoid is spherically 
aligned when (3 = 2 (in fact, this is a special case of Edding- 
ton's (1915) theorem). It is interesting that this somewhat 
exceptional case seems to reinforce our arguments for the 
physical nature of the spherically aligned Jeans solutions. 



6 6 CONCLUSIONS 

One way forward in the construction of three-integral galaxy 
models is to exploit partial integrals. These are invariants 
of the motion specific to particular orbital families. Phase 



space is therefore broken up into patches, and each patch is 
covered by an appropriate integral. This general picture is 
similar to that advocated by Binney and co-workers (e.g., 
Binney 1994; Kaasalainen & Binney 1994), who fit phase 
space tori to different orbital families using different map- 
pings. It also has some tradition in the theory of resonant 
orbits (Gerhard & Saha 1991; Dehnen & Gerhard 1993). The 
main result of this paper is the identification of such a partial 
integral for the thin and near-thin tubes in arbitrary axisym- 
metric scale-free potentials. It is a generalisation of the total 
angular momentum. By careful numerical orbit integration, 
the excellence of the partial integral for the thin tubes in one 
family of models - the scale-free power-law models (E94) - 
has been verified in detail. This supplements earlier work 
on the scale- free logarithmic potential (ZES). For all the or- 
bits we have examined in flattened potentials, this partial 
integral is always better conserved than the total angular 
momentum 

Approximate and partial integrals are needed in stel- 
lar dynamics for the construction of distribution functions 
(DFs). For the particular case of the scale-free power-law 
models, we carried out a detailed analysis of how to build 
three-integral DFs. This yielded the following conclusions: 

(1) If DFs depend on any approximate or partial integral, it 
is not guaranteed that they satisfy the Jeans equations. To 
be sure of building a galaxy model in which the motions of 
the stars balance the force field, it is important to combine 
both the Jeans and Boltzmann approaches. We have shown, 
both here and in ZES, that such combined Jeans and Boltz- 
mann solutions can be found in a simple way. 

(2) The solution of the Jeans equations for arbitrary ax- 
isymmetric scale-free models has been deduced. Jeans solu- 
tions in which the velocity ellipsoid is spherically aligned 
probably are physical, at least for moderate flattenings 
and anisotropies. This holds because the partial integral 
is quadratic in the velocities and does not generate any 
cross-terms in the stress tensor. At present, it is not known 
whether any of the remaining Jeans solutions are physical. 
Even at the level of the Jeans equations, it is evident that 
we cannot construct very radially anisotropic models with 
strong cusps {(3 > 0). 

(3) DFs depending on the classical integrals E and L z and 
the partial integral for near-thin tubes 73 can generate both 
the density and the spherically-aligned second moments of 
the power-law models. The scale-free power-law models are 
now a simple and useful set of galaxy models which have 
both two-integral and three-integral DFs readily available. 

There are three pressing questions left unresolved by the 
work in ZES and this paper. First, what is the appropriate 
partial integral to choose for the radial orbits in scale-free 
models? Although the phase space areas occupied by the 
thin and near-thin tubes are patched by the partial inte- 
gral deduced here, the areas occupied by the eccentric ra- 
dial orbits are not. The answer to this question is needed 
to build radially anisotropic three-integral DFs. Second, al- 
though one sub-set of the Jeans solutions displayed in Figs. 1 
and 2 has been shown to be physical by construction of suit- 
able DFs, the status of the remaining solutions is unclear. 
Is it possible that at least some - perhaps even all - of the 
remaining Jeans solutions are physically realisable? This is, 
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of course, related to the first question, as the Jeans solu- 
tions must be built by DFs depending on partial integrals 
appropriate to other orbital families. Third, how does this 
work generalise to fully triaxial scale-free models? Here, the 
solutions of the Jeans equations are now explicitly available 
(Carollo, de Zeeuw & Evans 1996), but there are no known 
integrals other than energy. 
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7 APPENDIX A: JEANS SOLUTIONS FOR 
CORED MODELS 

Here, we give the extension of solution (2.20) to the stresses 
associated with slightly more general axisymmetric densities 
in the potentials of the cored power-law models (E94). 
We consider the general density law: 



P = 



aiiT + az 2 + diR 2 c 



(R 2 + R 2 + z 2 /q 2 y+->/ 2 ' 
with oi, ci, d\ and 7 free parameters. The potential is 
< \\n{R 2 c +R 2 + z 2 /q 2 ), 

Rc 



~ P{Rl + R 2 + z 2 /q 2 f /2 ' 



for /3 = 0, 
for f3=£0. 



(Al) 



(A2) 



Here, R c is the core radius. When R c — we recover the 
scale-free potentials (2.1) with g(8) given in (2.4). The two- 
parameter solution of the Jeans equations is: 



p{vr) 



p(vrv z ) 



P( V z) 



p{vl) 



1 A!R 2 + C lZ 2 + D1R 2 
(a-1) (R 2 + z 2 /q 2 + R 2 C )<* ' 

1 B 2 Rz 

(a~T) (R 2 + z 2 /q 2 + R 2 ) a ' 

1 A A R 2 + C 3 z 2 + D3R 2 
(a-1) (R 2 + z 2 /q 2 + R 2 )<* ' 

_L A4R 2 + C1Z 2 + D!R 2 

(a-1) (R 2 + z 2 /q 2 + Rl)<* ' 



(A3) 



Here, a = (2 + 7 + (3)/2 and the relationships between the 
coefficients and the Hi and H2 parameters are given in Ta- 
ble 4. 

Self-consistent models have 

ai = ±(l-f3q 2 ), 
q 2 



ci 



di 



> 



1+/3n 



(A4) 



+ 2, 



7 = a = 2 + (3. 

Upon substitution of these expressions, and taking the limit 
R c = 0, we recover the results in Section 2.4 of the main 
text. 

It is instructive to transform to spheroidal coordinates 
(A, <f>,v), with (A, v) the two roots for r of R 2 /(t — R 2 ) + 
z 2 /(T-q 2 R 2 ) = 1. The foci lie &t(R = 0,z = ±R c y/l - q 2 ), 
and A and v are given by 



X,u= \[R 2 + z 2 + Rl{l + q 2 )]± 

\^{R 2 + z 2 ) 2 + 2R 2 {l-q 2 ){R 2 -z 2 ) + Ri{l-q 2 ) 2 . 



(A5) 
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Table 4. Coefficients for the Jeans solution (A3). 





[(a 


-l) + (2-a)q 2 ]ifi + (a-l)H 2 
. (a — l)ai +<j ci 






1 2a 


Ci 


(2- 




D\ 


r/ 

[(a 


-l) + (2-a)q z \H 1 


B 2 


(a- 


-i)^i 




(2- 


a)q 2 H 1+ (a 1 ^ +9 Cl 


c 3 


#1 




A 4 


(2- 


a)[2a-2+(3-2a)g 2 ]_f/i 

+(a 1)(3 2a)g 2 + 3 ^~ 1 ) ai + ( 3 ~ 2 ^' ?2ci 



The potential then simplifies to 

±lnA + \\nv + hiqR c , for/3 = 0, 



for 13 ± 0. 



(A6) 



The Jeans solutions (A3) can be rewritten to give p{v 2 }, 
p(v\v v ), and p(v 2 ) by means of the relations between 
(vr,v z ) and (v\,v u ) given in, e.g., Dejonghe & de Zeeuw 
(1988). In particular, we find 

H 2 Rz 



(A - v)p(v\v v ) = - 



{Rl + R 2 + z 2 /q 2 ) a ~ 1 ' 



(A7) 



This allows us to consider three limiting cases at once. 

When q = 1 the potential - but not necessarily the 
density - is spherical. The foci of the spheroidal coordinates 
now coincide with the origin, and hence (A, <f>, v) reduce to 
ordinary spherical coordinates, with v\ = v r and v v = — v$. 
It follows that (v r ve) oc H 2 - But in a spherical potential this 
cross term in the stress tensor must vanish by the symme- 
tries of the individual stellar orbits, so that the only Jeans 
solutions with physical DFs are those with H 2 = 0. This 
result is valid for arbitrary f3 and R c . 

When /3 = 2 the potential is of Stackel form in the 
coordinates (\,4>,u), for arbitrary q and R c (de Zeeuw & 
Lynden-Bell 1985). It then admits three exact integrals of 
motion that are quadratic in the velocity components, and 
any DF will therefore give (v\v„) = 0. In this limit we must 
therefore also restrict the Jeans solutions (A3) to those with 
H 2 = 0, irrespective of the values of q and R c . 

When R c = the spheroidal coordinates again reduce 
to sphericals, but the potential is flattened unless q = 1. In 
this case the model potential is scale-free, but there is no 
restriction on H 2 , unless q — 1 or j3 — 2. 

The above shows that if we consider our Jeans solutions 
as functions of the parameters q and (3, then along the two 
boundaries q = 1 and /3 = 2 of the parameter-space they 
must reduce to a one-parameter family, with H 2 = 0. There 
appears to be no physical reason for a similar restriction in 
the entire (q, /?)-plane. 

Finally, we remark that the line-of-sight projected sec- 
ond moments associated with the solutions (A3) are all of 
the general form (3.7) of Evans & de Zeeuw (1994), where 



the coefficients , bij , Cij and dij can be worked out explic- 
itly. 

8 APPENDIX B: CONSTRAINTS ON THE 
POSITIVITY OF THE PRINCIPAL STRESSES 

In this appendix, we write out explicitly the constraint that 
the two-parameter Jeans solutions discussed in Sections 2.4 
and 2.5 have positive principal components. The areas inside 
the following four curves must be deleted from Figs. 1 and 2. 

(l-Q)(Q 2 -Q+f3(Q 2 + Q+2H 1 )) 2 



2Q(l+f3)(2QH 1 -Q(Q-l) 2 +f3(2H 1 +Q-Q3))' 



„ (l-QKQ 2 -2Q-H 1 +/3Q 2 +/3H 1 ) 
H2 = q(T^ >< 

Q 2 ~Q+2(3H 1 +PQ(l+Q) 

Q{5Q-2Q 2 +4H 1 -3)-2H 1 +[3(2H 1 +Q+Q 2 -2Q3)' 

„ (l-Q)(Q+#iQ+/?(#iQ-2ffi-Q) 

H2 = WW) x 

f3{Q 2 +Q+2H 1 )-Q(l-Q) 

Q(2-2H 1 ~3Q+Q 2 )+f3(Q 3 +Q 2 +2H 1 Q-2Q-4:H 1 )' 

„ _ gi(i-Q) 

Q(l+/3)(2//!+g-Q 2 ) 2 
(Q(Q a -2Q 2 -Q~2H 1 (l+Q)+f3(Q 2 +Q 4 +4H 2 +4QH 1 )) 

±2Q[(2H 1 +Q)(Q-2PH 1 ~I3Q){2Hi+2Q-Q 2 ~PQ 2 )] 

These lengthy expressions are readily evaluated numerically. 
A quick check on their accuracy is provided by taking the 
limit (3 = 0, when they reduce to equations (2.38) of ZES. 



9 APPENDIX C: TRIAXIAL COMPONENTS 
WITH TRIAXIAL KINEMATICS 

Almost certainly, the potentials of galaxy haloes are not ex- 
actly spherical. But, one of the very attractive reasons for 
using spherical potentials is the ready availability of four 
functionally-independent integrals of the motion E,L x ,L y 
and L z . As realised clearly by White (1985) and Arnold 
(1992), this is useful for building axisymmetric tracer pop- 
ulations with triaxial kinematics - as is warranted by the 
shape of the local velocity ellipsoid in the Milky Way, for 
example. The same methods can also be exploited to give tri- 
axial tracer populations with triaxial kinematics! (c.f., Math- 
ieu, Dejonghe & Hui 1996). 

Together with White (1985), the work in Section 3 
shows that the DFs 



f m ,n(E, L 2 ,L 2 ) — 



l L 2n \E\ a , 



if p + 0, 



Vm , n L 2m L 2 z "e X p{-E), if = 0, 
generate the density law 

— *y • 2 rt n 

p — r 'sin 6, 

in the spherical scale-free power-law potentials 
llogr, if (3 = 0. 



(Cl) 



(C2) 



(C3) 
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The constant r\ m , n is available as (3.7) when (3 ^ 0. For the 
logarithmic case ((3 = 0), ZES show that 



'7,7 



( m+n+ 2)"+"+!?r(n+l) 
7iT(n+i)r(n+m+l) 



(C4) 



But, there is nothing special about the 2-axis, and no 
reason for the z-component of angular momentum to play 
a distinguished role! It is straightforward to show that the 
DFs 



g-m,n(E, L , L x ) — 



if p ^ 0, 



J?m , I1 X 2m L 2 ™exp(-.E), if /? = 0, 
give the stellar density 



(C5) 



p — r 



sin 2 9 sin 2 6 + cos 2 ( 



(C6) 



and the kinematics 



(vl) 



(ve) 



Voir 



2m + 2n + P + 7 



m + n + 1 (2n + 1) sin 2 <f> + cos 2 # cos 2 



n + 1 



sin 2 + cos 2 # cos 2 i 



'/ 2 \ 
-{v r ), 



1 21 m + w + 1 sin 2 4> + (2n + 1) cos 2 8 cos 2 </> , 2 i 



(C7) 



n + 1 



' + cos 2 # cos 2 



. n(m + ra+l) sin cos 6 cos 6 1 2( 
(vev*> = ——. . 2 7 , 1 „, a o , W-)- 



n + 1 
Equally, the DFs 



sin + cos 2 6> cos 2 <j> 



T Am r An I t-ii a 

/1 (V T 2 T 2 \ — i ? ? m >™ ij y 1^1 

ijm,n 



L' m L y n exp{-E), if/3: 



if P 0, 
0, 



(C8) 



correspond to 



sin 2 # cos 2 4> + cos 2 1 



(C9) 



with the kinematics (C7) on making the modifications 
cos <f> — * sin <f> and sin 4> — -> — cos <f>. Using the method of su- 
perposition of components, a very general DF corresponding 
to a triaxial halo with a triaxial velocity ellipsoid is 



f(E, L 2 ,L 2 , L 2 



m.n 



m.n 



(C10) 



The components can find a ready application for fitting the 
data on the kinematics of the stellar halo of the Milky Way. 
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